NASA TECHNICAL 

MEMORANDUM 


oo 


CM 

LT\ 


CO 


N 


66-13397 


(ACCESSION NUMBER) 


(NASA CR OR TMX OR AD NUMBER) 



L 


J2=- 


(CATEGORY) 


GPO PRICE $ 
CFSTI PRICE(S) $ 

Hard copy (HC) 
Microfiche (MF). 


'‘0V 




ff 653 July 65 


SATISFACTION OF ASYMPTOTIC BOUNDARY CONDITIONS IN THE 
NUMERICAL SOLUTION OF BOUNDARY-LAYER EQUATIONS 


by Philip R. Nachtsheim and Paul Swigert 
Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PREPRINT prepared for Ninth Midwestern Mechanics Conference, 
Madison, Wisconsin, August 16-18, 1965 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION - WASHINGTON, D.C. - 1965 


SATISFACTION OF ASYMPTOTIC BOUNDARY CONDITIONS IN THE 
NUMERICAL SOLUTION OF BOUNDARY -LAYER EQUATIONS 


by Philip R. Nachtsheim and Paul Swigert 

Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PREPRINT prepared for 
Ninth Midwestern Mechanics Conference 
Madison, Wisconsin, August 16-18, 1965 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


E-2959 


SATISFACTION OF ASYMPTOTIC BOUNDARY CONDITIONS IN THE 
NUMERICAL SOLUTION OF BOUNDARY-LAYER EQUATIONS 

-ft -ft-ft 

"by Philip R. Nachtsheim and Paul Swigert 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 


ABSTRACT 




A method for the numerical solution of differential equations of the 
boundary- layer type is presented. The asymptotic boundary conditions axe 
satisfied at the edge of the boundary layer by adjusting the initial con- 
ditions so that the mean square error between the computed variables and 
the asymptotic values is minimized. A least squares convergence criterion 
is used to propose a method of solution in which the edge of the boundary 
layer is approached in steps. The convergence rate to a solution is 
rapid and appears to be insensitive to the initial guesses of the initial 
conditions. Use of a least -squares convergence criterion leads to the 
unique solution even in the case of retarded main stream flows. A description 
of the method is given, and two examples taken from boundary- layer theory are 
worked . V 

INTRODUCTION 

The numerical integration of the ordinary differential equations of 
boundary-layer theory involves the ' satisfaction of asymptotic 


boundary conditions; that is, some boundary conditions are specified at the 
initial point or wall, and others are specified as limits that must be ap- 
proached at large values of the independent variable corresponding to the 
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edge of the boundary layer. In order to integrate the differential equations 
numerically from the wall to the edge of the boundary layer, it is necessary 
to specify as .many additional conditions at the wall as there are conditions 
to be satisfied at the edge of the boundary layer. These additional initial 
conditions have to be varied in order to satisfy the conditions at the edge 
of the boundary layer. All methods that have been used to find the proper 
initial conditions rely on the fact that, for large values of the independent 
variable, the integrals of the differential equations depend on the initial 
conditions. One method that has been used to find the proper initial con- 
ditions, called the initial -value method (ref. 1), is that of obtaining 
integrals of the differential equations with guessed initial conditions; 
that is, an attempt is made to integrate the differential equations to the 
edge of the boundary layer. If this can be accomplished, corrections are 
made for the initial guesses by the Newton-Raphson method, and the process 
is repeated until convergence is achieved. 

In order to carry out this iteration, additional differential equations 
have to be integrated. There will be as .many additional differential equa- 
tions as there are initial conditions to be varied. The integrals of these 
additional equations called perturbation equations give the rate of change 
of the integrals of the original differential equation with respect to the 
initial conditions. 

Another method that has been used, called quasilinearization (ref. 2), 
is similar in principle to the initial -value method, except that the dif- 
ferential equations are linearized. The resulting linear differential equa- 
tions for the current approximation are inhomogeneous, since they also contain 
the previous approximation as .members. 
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The perturbation equations of the initial -value method correspond to a 
complementary function of the linear differential equations of the quasi- 
linearization method. In application, the two methods are different in that, 
in the quasilinearization method, the complementary functions are used to 
obtain solutions to the inhomogeneous linear differential equations, whereas 
in the initial-value method, the solutions to the perturbation equations are 
used to adjust the starting values at the initial point. A disadvantage 
associated with the quasilinearization method is that functions representing 
the previous approximation have to be stored in order to construct the current 
approximation. Furthermore, the solution of the inhomogeneous equation ob- 
tained by combining complementary solutions and a particular integral is 
usually not well determined except near the origin or initial point, as pointed 
out by Hartree (ref. 3). The solutions of the perturbation equations, however, 
can be used to determine the proper initial conditions closely. This is the 
procedure that is followed in the initial -value method. 

Both methods have been used to obtain solutions of the boundary-layer 
equations. Neither of these methods, however, has solved the problem of when 
to stop the integration. Since cne boundary condition is specified as a limit 
that must be approached at large values of the independent variable, the inte- 
gration should be stopped at a value of the independent variable that is suffi- 
ciently large so that the various functions are approaching their asymptotic 
values. This value of the independent variable could properly be called the 
edge of the boundary layer. In the various methods used to solve this problem, 
this value of the independent variable must be guessed at beforehand. If the 
guessed value is too small, the integrals of the differential equation may not 
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"be able to satisfy the imposed conditions; or , there may be more than one 
value of an initial condition that leads to integrals of the differential 
equations that satisfy the imposed boundary conditions. 

If the guessed value is too large, it is possible that the integrals 
of the equations will diverge; or, that convergence to a solution is ex- 
tremely slow. 

Another problem that sometimes besets the numerical integration of the 
boundary-layer equations is the apparent insensitivity of the integrals of 
the boundary- layer equations to the initial conditions. This difficulty 
appears in the integration of the Falkner-Skan differential equation (ref. 4). 
All the integrals of the Falkner-Skan equation for the case of nonretarded 
flows in the main stream tend to diverge except the ones with the proper 
initial conditions; however, for the case of retarded flows, a solution with 
any value of the initial condition near the correct one will eventually meet 
the proper boundary condition at infinity. The unique solution was found by 
Hartree (ref. 4) for retarded flows by imposing the additional condition that 
the desired solution is the one that approaches the boundary condition most 
rapidly from below. However, this behavior of the solutions can only be 
learned from several trial runs. 

In an effort to adapt the initial-value method to the solution of dif- 
ferential equations with asymptotic boundary conditions, a method was developed 
to eliminate the problem of when to stop the integration. This method of 
solution is capable of satisfying the boundary conditions at the edge of the 
boundary layer correctly; that is, the boundary values are approached asymp- 
totically. This is accomplished by choosing the additional initial conditions 
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so that the mean square error between the computed variables and the 
asymptotic values is minimized. A least-squares convergence criterion is 
used to propose a method of solution in which the edge of the boundary 
layer is approached in steps. The convergence rate to a solution is rapid, 
and convergence appears to be insensitive to the initial guesses of the 
initial conditions. Use of a least-squares convergence criterion leads to 
the unique solution even in the case of retarded main stream flows. A 
description of the method is given and two examples taken from boundary-layer 
theory are worked. 

DESCRIPTION OF METHOD 

A description of the method can best be achieved by referring to a 
definite example of a boundary-value problem with an asymptotic boundary 
condition. In Falkner and Skan's treatment of the laminar boundary layer 
of an incompressible fluid, the following equation arises (see ref. 4) 

f" = - ff" + 0(f' 2 - 1) (1) 

with the boundary conditions 

T) = 0: f = f ' = 0 

T] -* oo{ f ' -» 1 

The primes denote differentiation with respect to q, a measure of the dis- 
tance from the wall. The dependent variable f is related to the usual 
stream function. The potential flow is given by a power law. The flow 
velocity outside the boundary layer in the main stream is proportional to 
the distance along the wall raised to the power 3/{2 - 3). For retarded 
main stream flow 3 < 0, and for nonretarded flow 3 > 0. 
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In practice, the asymptotic boundary condition is replaced by the 
condition that f ' = 1 to a sufficient degree of accuracy at r| = T| e £g e > 
where r^ge is the value of the independent variable at the edge of the 
boundary layer. The boundary-value problem is equivalent to the problem 
of finding a value of f"(0) for which the boundary condition at the edge 
of the boundary layer is satisfied; that is, it is desired to find a 
solution of the nonlinear equation 

[f"(0)J = 1 


f ' 

edge 


at H = 'Hedge* where f edge = f, (n e dge>* ^ function f edge[ f "( 0) ] 
in this problem is not in general explicit, it will be expressed as a 

function of f"(0) through an integration of equation (l). With the 
notation x= f"(0), observe that a small change Ax in x changes f' 
by the amount 

df 1 


5F- 


= Ax 


so that the necessary correction to a first approximation comes from the 
solution of the linear equation 


1 = f ' 


+ 



dx 


at H = Hedge* 

The solution of the equation for Ax can be performed provided the 
partial derivative of f ' with respect to x can be evaluated at "Hedge* 

The partial derivative can be evaluated by forming the perturbation equation 
for the function f'. This equation is obtained by differentiating the terms 
in equation (l) appropriately. With the notation 



f ' 
x 


df ' 


~ "c)x * 


etc. 
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The following perturbation equation is obtained 

f^" = - (ff£ + f"f x ) + 2pf'f^ (2) 

with the initial conditions 

’l - 0; f X - f i - °> f x > 1 

Given an initial estimate of x = f"(0), the subsequent values of x 
are computed automatically in the iteration procedure described previously. 

The problem of where to stop the integration remains; that is, q e £ ge 
is unknown. Figure l(a) illustrates the difficulties that arise when the 
boundary condition f ' = 1 is applied at a finite value of q. 

The data for figure l(a) were obtained from calculations that are 
described in detail subsequently. For the purpose of understanding the 
results presented in this figure, the reader should know that equation (l) 
was integrated for specified values of f"(0) with 2=1. The integration 
was stopped at q = 5. Tentative values of q edge will hereinafter be 
referred to as q s top‘ 

Figure 1(a) shows f' and f " as a function of x = f"(0) at 
q s top = It can be seen that there are two values of x that satisfy 
the condition f ' = 1 when il g ^ 0 p I s taken to be 5, namely x = 0.85 and 
x = 1.23. Figures 1(b) and (c) show the two solutions corresponding to the 
two values of x. The proper value of x can be determined by observing 
where f" = 0. From figure 1(a) it can be seen that the correct value of 
x is approximately 1.23. » It can be concluded from this illustration that 
the satisfaction of the boundary condition f’ = 1 at q s -top = 5 does not 
ensure the satisfaction of the asymptotic boundary condition on f', namely 


8 


that the value of f ' approaches the value 1 asymptotically. In order to 
satisfy the asymptotic condition at a finite value of r) , the condition 
f 1 = 1 and f" = 0 should both he imposed (examination of the differential 
equation (l) will reveal that all higher derivatives will be zero also); 
that is, Ax should be chosen so that both equations 


1 = f ' + f i Ax 


0 = f" + Ax 


(3) 

(4) 


are satisfied at t| = Tig-top* This is in general impossible since there are 
two conditions and only one adjustable parameter, Ax. However, a satisfactory 
solution that is consistent with the idea that the boundary condition cannot be 
satisfied exactly at a finite value of r\ is to seek the least -squares solution 
of the preceding equations (ref. 5); that is 

fi d - f ') - f,,f x 


Ax = 


(5) 


+ f ■ 


x x 

In this way all the information contained in equations (3) and (4) is 
retained. 

The least-squares solution of equations (3) and (4) attempts to minimize 
the sum of the squares of the errors in these equations as a function of Ax. 
The value of x that gives Ax = 0 corresponds to the minimum with respect to 
x of the quantity 


E = (1 - f') 2 + f" 2 (6) 

as can be seen from equation (5). In figure 2, the quantity E is shown 
plotted as a function of x for various values of T)stop* The ^ or 

figure 2 were also obtained from the calculations described subsequently. 
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For this example (p = 1 and "n s -t 0 p = 5) it is seen that the least -squares 
solution gives a unique solution for x. It is this property, that the 
minimpm of the sum of the squares of the appropriate boundary functions 
corresponds to no change in the required initial conditions, that enables 
satisfaction approximately of the asymptotic boundary condition at a finite 
value of the independent variable. The edge of the boundary layer could 
be defined as that value of r| for which the minimum E is less than some 
preassigned small value as the range of integration is increased. This gives 
a value of q edge . 

The concept of convergence in the least-squares sense permits modification 
of the initial-value method in order to solve a wide class of boundary-layer 
problems. To understand why the initial value method requires modification, 
consider how the method can fail when it is applied to boundary-layer type 
problems . 

Recall that in the initial -value method, corrections are obtained to the 
guessed initial conditions after integrating the differential equations to 
the edge of the boundary layer. This method can fail if the initial guess is 
so poor that the integration "blows up" (i.e., diverges beyond a prescribed 
limit) before reaching the edge of the boundary layer. 

This problem can be avoided and corrections to the initial conditions 
can be obtained by attempting to minimize the mean square error between the 
computed solution and the asymptotic values before reaching the edge of the 
boundary layer. The modification of the initial -value method consists then 
of carrying out the minimizing process at any arbitrary value of tig-top* When 
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this modified procedure is used, the choice of ri s ^ 0 p^ where the minimizing 
process is first carried out, may he dictated by other considerations than 
where the edge of the boundary layer is, such as keeping the variables 
within prescribed limits so that the solution does not diverge. After the 
.minimum has been found and corrections to the initial conditions made at 
the chosen value of Ti s ^ 0 p, the range of integration can be increased if 
necessary. 

If the value of the minimum found is not small enough, it will be 

necessary to increase the range of integration. In this way, the edge 

of the boundary layer is approached in steps by successively increasing 

the range of integration until the correct initial conditions and the value 

of n , are found. 

'edge 

For this example, the Falkner-Skan solution with 3=1, it appears 
the quantity E has a definite minimum at each value of q including 
T] = 0 as can be seen from figure 2. The numerical solution of this example 
will be described in the next section. 

NUMERICAL SOLUTION 

The method described in the previous section was programed in 
FORTRAN IV, double precision on the IBM 709411-7040 direct-couple system. 

The boundary-layer and perturbation differential equations were rewritten 
as systems of first-order differential equations and integrated with a 
predictor-corrector (Adams -Moulton) subroutine using one correction per 
step and a fixed increment. 

The boundary-layer equations along with the perturbation equations 
were integrated to the specified value of the independent variable. The 
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corrections were then determined, and the process was repeated until the 
relative change in the correction term equation (5) was less than a small 
preassigned value. 

Extreme accuracy was not required at the smaller values of xi s -top 
since the boundary conditions cannot be well satisfied there. A value of 
1x10’® was used at the larger values of ii s -top check the relative 
change in the correction term. This small value seemed reasonable since 
convergence is so rapid. 

RESULTS OF NUMERICAL SOLUTION FOR EXAMPLE 1 

Figure 3 shows the solution f’ of equation (l) with 3=1 as 
a function of the independent variable t]> for the first five trials needed 
for convergence. With an initial guess of 1.0 for f"(0), the equations 
were integrated to a value of rj s ^ 0 p = 2 twice and to a value of 
Tstop = 5 ' t * iree times. Observe that, for the first integration, the 
values of f 1 tend to deviate radically from the correct solution at a 
relatively small value of t). If these calculations had been carried out 
to a larger value of r\, the values of f ' would have attained a magnitude 
so large that any use of these large numbers in a Newton-Raphson scheme would 
be meaningless. 

The effect of ri s ^ 0 p on the convergence process is illustrated in 
figures 2 and 4. Figure 2 shows the error, equation (6), as a function of 
f"(0) for different values of ii s top* Observe that by initially using a 
small value for (r| sto p = 2 in fig. 2) the range of meaningful initial 

guesses is increased, that is, initial guesses that yield a relatively small 
error. For extremely large values of T] s ^ 0 p 


the parabola-like curves in 



12 


figure 2 degenerate to a vertical line, and the initial guess is limited to 
the correct answer. 

Figure 4 shows the effect of ri s ^ 0 p on the convergence rate. The value 
of f"(0) after two trials is plotted against the initial guess of f"(0) 
for different values of Hstop" ® ie distance of a point on these curves from 
the straight line shown in the figure is a measure of the rate of convergence. 

If a point remains on this line after two trials, there is no tendency toward 
convergence at all. It can he seen from figure 4 that the points on the 
curve for %-fcop = 2 are a l wa Y s farther from the straight line than the points 
on the curve for Ti s ^ 0 p = 5. This indicates that the rate of convergence is 
.much greater for T^gtop = 2 than for = 5 for all initial guesses of 

f"(0). 

To give an idea of the rate of convergence to a solution for this example 
using an initial guess of f"(0) = 1, table 1(a) is presented. Using quasi- 
linearization and five corrections Radbill (ref. 2) obtains a value of f"(0) 
accurate to one decimal place compared with Yohner and Hansen (ref. 6), whereas 
the .method developed herein reached an answer after two corrections that agrees 
with that of Yohner and Hanson to two decimal places. After two more corrections, 
the value of f"(o) agrees with that of Yohner and Hansen (ref. 6) to within 
2 units in the seventh decimal place. The last integration was performed as a 
stopping criteria for the computer. 

Computer time for the complete solution, using a step size of 2 -4 was 
approximately 0.05 minute. 

This example was also programed by single precision with good results. 

The difficulty in using single precision presented itself in the form of reducing 


13 


the range of initial guesses that could successfully be used to obtain the 
correct answer. 

Equation (l) was also integrated for the case of a retarded main stream 
flow with |3=-0.1. However, no additional conditions had to be imposed 
such as those used by Hartree (ref. 4) in order to obtain the unique solu- 
tion. The value of f"(0), which was obtained by minimizing. the mean 
square error, agreed with the value obtained by Smith (ref. 7) for this 
case to six decimal places. 

EXTENSION TO SYSTEMS OF EQUATIONS 

The least squares modification of the initial-value method described 
previously for a differential equation with a single dependent variable can 
easily be generalized to systems of equations. As an example of a problem 
with two dependent variables, consider the boundary value problem for the 
free -convection flow about a vertical plate. This example is taken from 
reference 8. It consists of solving the differential equations 


f"’= 3ff" + 2f' 2 - 

h 

(7) 

h" = - 3Prfh' 


(8) 

with the boundary conditions 



q = 0: f = f ' = 0, 

h = 1 


q ^ oo: f -* 0, h 

-* 0 



Again the primes denote differentiation with respect to q (a measure 
of the distance from the wall), and the dependent variable f is related 
to the usual stream function. The dependent variable h is proportional 
to the temperature excess of the fluid over the ambient temperature, and 
Pr denotes the Prandl number. 
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f ' 
edge 

h, 


Since there are two asymptotic boundary conditions to satisfy, two 
additional initial conditions at the wall will have to be adjusted; that 
is, values of x = f"{0) and y b h'(0) are sought that will satisfy 
simultaneously the nonlinear equations 

f " (0) ,h ' (0)j = 0 

edg e [ f, '(0);h , (0)j =0 

at n = T} , where f’ = f'(ri ) and h = h(r) ). The 
edge edge edge edge edge 

necessary first corrections to a first approximation x and y come from 

a solution of the linear equations 

0 = f* + f' Ax + f' Ay 
x y ^ 

0=h+h x Ax+hyAy 

at t) = q ^ . Where f x = , f ' = etc. However, in order to 

B Sx ^ dy 

satisfy the asymptotic boundary conditions, the preceding equation must be 
supplemented by 

0 = f" + f x Ax + f" Ay 
0 = h' +h x Ax+hyAy 

and Ax and Ay must be found such that the sum of squares of the errors 
in the preceding four equations be a minimum. The least-squares solution 
for those quantities is given by the solution of the following matrix equation: 



f '2 + h2 + f "2 + h ' 2 
XXX X * 

f'f' + h h + f"f" + h'h' 
x y x y x y x y 


Zbc 

f 'f ' 
X y 

b. 

+ h h + f "f " + h'h' , 
x y x y x y’ 

f . 2 + h 2 + f + h . 2 

r 


Ay 

- 


f 'f ' + hh + f"f" + h'h' 

-X- X X X 

f 'fy + hhy + f"fy + h'hy 


( 9 ) 


The values of x and y that give Ax = Ay”‘ ; * 0 correspond to the 

* 

minimum with respect to x and y of the quantity 
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E = f ' 2 + h 2 + f" 2 + h' 2 (10) 

The partial derivatives with respect to x and y that appear in 
the equation (9) are obtained by integrating the appropriate perturbation 
differential equations. The perturbation differential equations for the 
x derivatives are 

f"'= - 3(f f" + ff") + 4f'f' - h (11) 

x x x xx 

and 

h x = - 3 Pr(f x h ' + fh x ) (12) 

with the initial conditions 

q = 0: f x = f x = h x = h x = 0 j f x - 1 

The perturbation differential equations for the y derivatives are 


= - 3(fyf" + ffj) + 4f 'fy - hy 

(13) 

hy = - 3Er(fyh' + fh£) 

(14) 


with the initial conditions 

n = 0: f y = fy = fy = h y = 0) h^. = 1 

Note that the equations for the y derivatives are the same as the 
equations for the x derivatives. The integrals of these equations will 
differ since the initial conditions are different. 

Note also that there are as many additional systems of perturbation 
equations to integrate as there are asymptotic boundary conditions to .meet. 

NUMERICAL SOLUTION OF EXAMPLE 2 AND RESULTS 
The same general procedure used in the previous example was employed 
to solve this example. One notable difference is that five trials were 
needed at the larger value of T]stop> rather than the three trials of the 
first example. This is probably because of the need of adjusting two initial 
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conditions rather than one, and because an error in one initial condition 
leads to an error in computing the correction term of the other. 

The solutions of equations (7) and (8) are shown in figure 5 for 
Fr = 0.733. Note (as in the previous example) that the functions f' and 
h of figure 5 diverge radically from the true solutions, for the initial 
guesses of trial 1. 

Table 1(b) shows the rapid rate of convergence for this more difficult 
example. With initial guesses of 1.0 and - 1.0 for f"(0) and h'(0), 
respectively, convergence is realized after two trials with = ^ an< ^ 

five trials with r| s -t 0 p = 8* The results are in close agreement with the 
published values of Ostrach (ref. 8). 

The effect of r| s -k 0 p on the convergence process for this example 
is illustrated in figures 6(a) and (b). Level curves for the error, equa- 
tion (10), are plotted in the f"(0) plane for r| s top = 2 in figure 6(a) 
and for t s top = 5 in figure 6(b). As in the previous example, by using 
a small value for T l shop ^ ran § e initial guesses that yields a 

relatively small error is again increased; that is, the area where E is 
less than 1 is much greater in figure 6(a) than in figure 6(b). For 
extremely large values of Ustop-’ the area enclosed by the level curve 
where E is less than 1 shrinks to a point, and the initial guess is limited 
to the correct answer. 

Computer time for the second example was about 0.15 minutes with a 
step size of 2”^. In both examples, the convergence was made insensitive 
to the initial guesses by choosing a small value of f° r the fi rs t 

two trials. The answer thus obtained for the initial conditions was good 
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enough to step Tig-top U P considerably. The error term may be .made as 
small as desirable by continuous stepping of Tig-top* This error term may 
also be used as a stopping criteria for computer calculations and findings 
Tiedge automatically. 

CONCLUSIONS 

The two main problems of integrating the boundary- layer equations, approx- 
imating the missing initial conditions, and finding "n e a.ge , * iave ^ een reduced 
to an automatic initial -value technique that is easily programed on high-speed 
computers. The method appears to be insensitive to initial guesses and con- 
verges quickly to the solution. 
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TABLE I. - CONVERGENCE HISTORY OF 
INITIAL VALUES 


(a) Example 1 


Trial 

Correction 


f"(0) 

1 


2.0 

1.0 

2 

1 

2.0 

1.2463981 

3 

2 

5.0 

1.2266764 

4 

3 

5.0 

1.2326729 

5 

4 

5.0 

1.2325878 


5 


1.2325878 


Published values 


Source 

Si 

f"(0) 

Radbill (ref. 2) 
Yohner and Hansen 
(ref. 6) 

5.0 : 

4.0 

1.2397 

1.2325876 


(b) Example 2 


Trial 

Correction 


f"(0) 

h'(0) 

1 


2.0 

1.0 

-1.0 

2 

1 

2.0 

.61735605 

-.59008297 

3 

2 

8.0 

.63211384 

-.57655550 

4 

3 

8,0 

.66795925 

-.52593128 

5 

4 

8.0 

.67387305 

-.50916398 

6 

5 

8.0 

.67412049 

-.50790411 

7 

6 

8.0 

.67412438 

-.50789273 


7 


.67412438 

-.50789273 


X ^ 


Published values 
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Figure 1. - Variation of f and f'. 






Figure 3. - Convergence history of solution of equation (1) 
for (5 * 1. 
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Figure 6. - Level curves of error as function of F'(0) and h'(0). 
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